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Theoretical and numerical studies are presented of the nonlinear amplitude modulation of dust- 

, acoustic (DA) waves propagating in an unmagnetized three component, weakly-coupled, fully ionized 

■ plasma consisting of electrons, positive ions and charged dust particles, considering perturbations 

' oblique to the carrier wave propagation direction. The stability analysis, based on a nonlinear 

fSJ , Schrodinger-type equation (NLSE), shows that the wave may become unstable; the stability criteria 

. depend on the angle 9 between the modulation and propagation directions. Explicit expressions for 

' the instability rate and threshold have been obtained in terms of the dispersion laws of the system. 

I The possibility and conditions for the existence of different types of localized excitations have also 

, been discussed. 
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PACS numbers: 52.27.Lw, 52.35.Fp, 52.35.Mw, 52.35.Sb 



I. INTRODUCTION 



The study of the dynamics of dust contaminated plasmas (DP) has recently received considerable interest due to 
their occurrence in real charged particle systems, e.g. in space and laboratory plasmas and the novel physics involved 
Q in their description 0. An issue of particular interest is the existence of special acoustic-like oscillatory modes, e.g. 
' the dust-acoustic waves (DAW) and dust-ion- acoustic waves (DIAW), which were theoretically predicted about a 
decade ago 0, Q and later experimentally confirmed 0, Q. The DAW, which we consider herein, relies on a new 
, physical mechanism in which inertial dust grains oscillate against a thermalized background of electrons and ions 
O i' which provide the necessary restoring force. The phase speed of the DAW is much smaller than the electron and ion 
thermal speeds, and the DAW frequency is below the dust plasma frequency. 
I ' A long-known generic characteristic of nonlinear wave propagation is amplitude modulation due to the nonlinear self- 
^ [ interaction of the carrier wave, which generates higher harmonics. The standard method for studying this mechanism 
is a multiple space and time scale technique I^Qj which leads to a nonlinear Schrodinger-type equation (NLSE) 
describing the evolution of the wave envelope. It has been shown that, under certain conditions, waves may develop a 
■ Benjamin- Feir- type (modulational) instability (MI), i.e. their modulated envelope may collapse under the influence of 
external perturbations. Furthermore, the NLSE-based analysis, already present in a wide variety of contexts |3.l^ll0|. 
reveals the possibility of the existence of localized excitations (solitary wave structures) whose form and behaviour 
depends on criteria similar to the ones necessary for the MI to occur. 

Not surprisingly, plasma wave theory has provided an excellent test bed for this approach since a long time ago 
[7, 11 - 19] and dusty plasma waves were no exception [20 - 22]. Among other noteworthy results, electron plasma 
modes have been shown to be stable against parallel modulation so do the ion plasma modes, yet only for 
perturbations below a specific wavenumber threshold [l^. Electron and ion acoustic modes, even though stable to 
r~ I ' parallel modulation 0, 0, are found to be unstable if one takes into account finite temperature effects 
O |. |12.ll5l ll q] or, most interesting to us, when subject to an oblique modulation of the wave amplitude [17 - 19]. These 
ILJ ■ results, based on Poisson - moment plasma equations, have been confirmed by similar studies from a kinetic point 
. ^ I of view |24| , for the ion - acoustic wave in an electron - ion plasma. In dusty plasma, the amplitude modulation of 
the DAW and DIAW has been investigated in Ref. [20 - 22]; similar studies have been carried out for oscillations 
in (strongly-coupled) dusty plasma quasi-crystals [25 - 26] . Finally, let us mention that attempts have been made to 
. . . refine the description of the DIAW modulation by including non-planar geometry effects ,27] , following an idea app lied 
earlier in the KdV (Korteweg-de Vries) description of a dusty plasma [23 an d dust-charge fluctuation effects [21|, an 
issue of particular importance in the present-time DP surveys (see e.g. |22,|33; also |l|). These effects, omitted in 
the present investigation, will be considered in a forthcoming work. 
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In this paper, we study the modulational instabihty of dust-acoustic plasma waves propagating in an unmagnetized 
plasma contaminated by a population of charged dust grains, whose dimensions and charge are assumed constant, 
for simplicity. Amplitude modulation is allowed to take place in an oblique direction, at an angle with respect to 
the carrier wave propagation direction. Once an explicit criterion for the occurrence of instability is established, our 
aim is to trace the influence of 6 on the conditions for the MI onset, and determine the magnitude of the associated 
instability growth rate. Finally, we shall also examine the possibility of the formation of localized excitations and 
discuss their characteristics. Exact new expressions are derived for all quantities of interest, in terms of the system's 
dispersion laws. Among other physical parameters discussed, our formulation leaves open the choice of sign of dust 
charge {signqd = ±1) (most often taken to be negative since this is the most frequently occurring case (l|) and the 
dust pressure ('temperature') scaling. Our aim in doing so is to address, among others, the question of the influence 
of the dust charge sign on the amplitude modulation mechanism. We may also attempt to clarify the effect of taking 
(or not) into account the dust pressure evolution equation (omitted e.g. in 20]) in the analysis. 

The manuscript is organized as follows. In the next Section, the analytical model is introduced. In Section III, 
we carry out a perturbative analysis by introducing appropriate slow space and time evolution scales, and derive a 
NLS-type equation which governs the (slow) amplitude evolution in time and space. The exact form of dispersion and 
nonlinearity coefficients in the NLS-type equation is presented and discussed. In Section IV, we carry out a stability 
analysis of the NLSE allowing for a thorough study of the DAW stability in various regions of the physical parameters 
involved. The analysis is pursued in Section V, where we discuss the possibility of the existence of localized solutions 
of the NLSE, and identify their forms in different parameter regions. Finally, we briefly summarize our results in the 
concluding Section. 



II. THE MODEL 



We consider a three component collisionless unmagnetized dusty plasma consisting of electrons (mass m, charge e), 
ions (mass m^, charge qi = +Zie) and heavy dust particulates (mass m^, charge qd — sZdc), henceforth denoted by 
e, i, d respectively. Dust mass and charge will be taken to be constant, for simplicity. Note that both negative and 
positive dust charge cases are considered, distinguished by the charge sign s = sgnqd = ±1 in the formulae below. 



A. Evolution equations 



The basis of our study includes the moment - Poisson system of equations for the dust particles and Boltzmann 
distributed electrons and ions. The dust (number) density nd is governed by the (continuity) equation 



^-fV.(n,u,) = 0, (1) 



and the dust mean velocity obeys 



^ + u,.Vu, = -*Lv<i>-^Vp,, (2) 

at rUd rridUd 

where $ is the electric potential. The dust pressure pd obeys 

+ Urf • Vpd = -7Pd V • Ud (3) 

Here 7 (/ + 2)// is the ratio of specific heats (/ is the number of degrees of freedom) e.g. 7 = 3 in the adiabatic 
one-dimensional (Id) case and 7 = 2 in the two-dimensional (2d) case. The system is closed with Poisson's equation 

V^* = -47r ^ = All (ji^e- qiUi- qdUd) = e{ne- ZiUi- s ZdUd); (4) 

note that the right-hand-side cancels at equilibrium due to the overall neutrality condition 

ne,o e - Uifl qi - Ud.o qd = ■ (5) 
The right-hand side in is often formulated in terms of the ratio fi = ne.a/{Zini,o): for convenience, we have 

fi = 1 + s- (6) 
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due to lO, so that a value lower (higher) than 1 corresponds to negative (positive) dust charge; fi obviously tends 
to unity in the absence of dust (in any case, /i > 0). We will retain this notation in the following, for the sake of 
reference to previous works. 

The electrons and ions are assumed to be close to a Maxwellian equilibrium. The corresponding densities are 

and 

n.«n,,oe-^'^*/'=-^', (7) 
where denotes the temperature of species s = e,i {ks is the Boltzmann constant). 



B. Reduced equations 

Re-scaling all variables over appropriately chosen quantities and developing around <i> = 0, Eqs. CJ - Q) can be 
cast in the reduced form 

dn 

-+V.(nu) . 0, 

— + u-Vu = -sV(/) Vp, 

at n 

dp ^ ^ 
— + u-\/p — — 7pV-u, 
at 

and 

V^^ = </)-a(/)2 + a'0^-s/3(n-l), (8) 

where all quantities are non-dimensional: n = Ud/ridfi, u = Ud/vQ, p = Pd/po and cj) = ^/^o; the scaling quantities 
(index 0) are, respectively: the equilibrium density nd,o, the 'dust sound speed' vq = Vd = {kBTe/md)^^'^ , po = 
n-dfiksTe and $o = {ksTe/Zde). Space and time in © are, respectively, scaled over: the DP effective Debye length 
'^D.eff = {^^e + '^Iji)^^^'^ (where Xd.s = (ksTs/'i'Knsfiqiy^^, s = e,i) and the inverse DP plasma frequency 
^pd ~ (47i'?^d,o'7d/W(i)~^/^. Recall that s =sgn qd, so the influence of the dust charge sign will be traced via the 
appearance of s in the forthcoming formulae. Finally, a — po/ {udflkBTe) is equal to unity, given the above choice for 
Po; nevertheless, a - often interpreted as a temperature ratio via a different scaling, see e.g. |2lj| - will be retained in 
order to 'tag' the influence of the coupling to pressure evolution equation (jSj) being taken into account - as compared 
to a previous work _20J where Eq. ^ has been omitted. As a matter of fact, expressions (9) - (11) therein are readily 
recovered here upon setting cr = 0, s = —1, a' = in Eq. (jHJl above. 
The dimensionless parameters appearing in ^ are 

1 ^» ^T-) ~ ~ , _ 1 ^» \t-) ~ + ^ 

^ T-~ + ^ d A + t 

and 

V Vd J \Vd 

where cd — ^Def f^p.d is the DA speed 0. Alternatively, in terms of /i defined above, one has: a — — ^ l^a^^+a^ ; ~ 
ijl^rrS' ^ ^ai+aa ' "^^cre Qfi = ^ and 02 = |^ ^ . AU these parameters are positive For /i < ^ = ^, 
we have the approximate expressions: a k, ^^ = 2^i^ ^^'^ a' ^ — ■§^'^ ~ | a^, as in also: (5 « 

A comment should be made, regarding the order of magnitude of the parameters a, a', (3. Notice that a takes very 
small (positive) values (as low as, say, 10~^ to 10~^) and so does a'; however, /3 may take high values, e.g. ranging 
from zero (for fJ. = 1 i.e. no dust) to, say, 10'^ — 10'^. Therefore, the numerical result of the scaling in our (DAW) case 
is completely different from the one in the dust ion-acoustic (DIAW) case despite the apparent similarity in the 
model expressions 113, H^; this is why we chose not to analyse the DIAW case any further, in the same text. 
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III. PERTURBATIVE ANALYSIS 
A. Outline of the method 

Let S be the state (column) vector (n, u, p, 0)"^, describing the system's state at a given position r and instant t. 
We shall consider small deviations from the equilibrium state = (1, 0, 1, 0)^ by taking 

OO 

S = + e ^(1) + 5(2) ^ ^ ^(0) + ^ , 

n=l 

where e ^ 1 is a smallness parameter. Following the standard multiple scale (reductive perturbation) technique 
we shall consider the following stretched (slow) space and time variables 

= e{x- Xt), T = e^t, (9) 

where A, bearing dimensions of velocity, is to be later interpreted as the group velocity in the x direction. In order to 
take into account the influence of an oblique amplitude modulation on the DA wave, we will assume that all perturbed 
states depend on the fast scales via the phase 9i = ]<. ■ r — cut only, while the slow scales enter the argument of the 

(n) 

l—th. harmonic amplitude , which is allowed to vary along x, 

oo 
/ — — oo 

The reality condition si""* — Si"''' is met by all state variables. Note that the (choice of) direction of the propagation 
remains arbitrary, yet modulation is allowed to take place in an oblique direction, characterized by a pitch angle 6. 
Assuming the modulation direction to define the x— axis, the wave- number vector k is taken to be k = [k^, ky) = 
{k COS0, k sin6'). According to these considerations, the derivative operators in the above equations are treated as 
follows 



d d d 2 ^ 

dt'^ dt~^ dC^^ d^' 



and 



i.e. explicitly 



V-> V + ex — , 



a2 a2 



and 



y2 , («) ^^W, ^ _ ;2^2 An) ^ 3^ ^^^^ ^ ^2 \ ^UO, 



for any A["'^ of the components of S'|"'' . 
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B. Amplitude evolution equations 

By substituting the above expressions into the system of equations ^ and isolating distinct orders in e, we obtain 
the nth-order reduced equations 



ilujn 



(») 



(n-l) 



dn 



(n-2) 



oo 

E E 



dr 



dc 



oo oo 



9u 



1 i' = -c 

(n-l) 
I 



(n-2) 



5C 



9r 

OO OO 

E E 

n' — 1 oc 



oo oo 

E E 



(n — n') 



a (n'-l) 

— liLOU), ' + S Zi k0,\ ' — X — 



ft ^^'^ 

77'k • + ,/"-"'-!) ^"i' 



+ C7 ( i/p|"^k 



(«'-2) 



9u 



oo oo 

E E 

n" — 1 oo 



dC dr 

(n'-n"-l) 



'■{l'-l"),x 



dC 



= 0, 



(10) 



(11) 



(n-l) 



(n-2) 



du 



(n-l) 



dr 



1 ■ 



dC 



oo oo 

E E P'-/ 

n' — l l' — — oo 



oo oo 



n' = l i'=-oo ^ ^ 



0, 



(12) 



and 



(Pfc^ + 1) + 



s /3 + 2ilk, 



(n-l) 



2j,(»-2) 



n' — l — — oo 



oo 



oo 



E E 



a' E E 

n',n" — 1 V ,1" —~oo 



(n—n' — n") An') /(n") 



= 0. 



(13) 



Notice the last three lines in Eq. Hll|) . which are due to the consideration of the pressure evolution equation and 
are absent e.g. in Ref. |23| - cf. Eq. (33) therein. Even though it is a which introduces coupling to p2|) (which 
becomes decoupled from the rest, that is, for a = 0), should one correctly consider the limit cr = in Eqs. ©, one 
should discard all three of the last lines in 111|) . For convenience, one may consider instead of the vectorial relation 
(|ll|l the one obtained by taking its scalar product with the wavenumber k. Finally, we see that Eqs. (32) - (34) of 
Ref. ^] are readily recovered upon setting cr = 0, s = — 1 and a' = in the above relations. 



C. First order in e: first harmonics and dispersion relation 

The first order {n = 2) equations read 



iluon 



(1) 



-ilujxi 



(1) 



ilk ■ u[^^ = , 



silk(l)[^^ 



(14) 
(15) 
(16) 
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and 

- (/2fc2 + 1)0(1) ^ 0. (17) 

For 1 = 1, these equations determine the first harmonics of the perturbation. The foUowing dispersion relation is 
obtained 

.^=^^+jak\ (18) 
Restoring dimensions, one may easily check that the standard DAW dispersion relation is thus exactly recovered: 



2 2 , fcsTd 2 _ fc2 2 7 2 /1 n^ 



The first harmonic amplitudes may now be expressed in terms of the first order potential correction 0^^^; we obtain 
the relations 

(1) ^ /) (1) l + fc^w (1) (11) ,(1) 

' ~ ■^*^°^^'^i ^ — ^3 — — cost^^i'^Cg (Pi , 

u^^l = i^sin0n« = .i±^i^sin^0W, (20) 

retaining, for later use, the (obvious) definitions of the coefficients c^-^^-* (j = 1,...,5) relating the state variables to 
the Ist-order potential correction (jy'p (so c'l^^ = 1). 

D. Second order in e: group velocity, 0th and 2nd harmonics 

The second order (n = 2) equations for the first harmonics provide the compatibility condition: A = Vg{k) = = 
Li;'(A;) COS0 = ^ [{T+I^p' + 7'''] cos^; the group velocity Vg can be cast in the form 

^3 ^ + ^^(l + fc2)2 ^3 

^ F [/3 + a7(l + F)]2 ^'^^ 



where we have denoted 



"^-'^[/3 + a7(l + P)]2- ^^^^ 



Note that vi ^ 1 in the limit cr 0, recovering exactly Eq. (43) in Ref. |20|. 
The 2nd-order corrections to the first harmonic amplitudes are now given by 

= ^.i[i(l + fc2)_2fccos0]^^zcr)^, 

1 (2) (2) , ;2n/^ ^ /)^ ^"^^l^' ■ (21) 5'?^1^^ 

(2) (2) _ . (21) d4>^i'^ 

Pi ^ = %c\ 



.(2) _ ,^5^ 



a(i) 



^1 



5C 



7 



and 



''il 



I s- 



7 



(21) d(j)[^^ 

C5 



cos 9 + [Vg—cos6-aj 



f3 



to 



/3 



9C 



9C 



(23) 



The choice of the value of A is arbitrary; we shall take A — Q. 

The equations for n = 2, / = 2 provide the amplitudes of the second order harmonics, which are found to be 
proportional to the square of the corresponding elements e.g. in terms of 0^^^ 



,(2) 



k • 



(2) 



(2) 
P2 



1 



-A 



(l + fc2)2 



J22) ,(1)^ 



A 



6/33fc2 



2sal3' + 3/3(l + fc^)(l + 2fc^)+ 27V (1 + /c")^ (1 + 4fc") 



1. + 



= ,(1)^ 



and 



if 



4fc2 



1 



s/3 



1a 



/?2 



, ,(1)2 _ (22) ,(1)2 

a\(j)\' = c\ ' (t)\' . 



(24) 



Notice that these expressions are isotropic i.e. independent of the value of 9. 

The nonlinear self-interaction of the carrier wave also results in the creation of a zeroth harmonic, in this order; its 
strength is analytically determined by taking into account the I = component of the three first third-order reduced 
equations (i.e. (|10|) - (|12l) for n = 3, / = 0) together with the corresponding fourth 2nd-order equation (i.e. (|13|l for 
n = 2, I — Q). The result is conveniently expressed in terms of the square modulus of the {n ~ \, I ^ \) quantities, 
e.g. in terms of ~ ("^(^H* ^(1) 



.(2) 



■'I )* VI 

-1 



1 



/3 + 7cr - ^2 j3 



1 + 2sal3 +k'^ +2 cos^ 
(1 + P)2 



-1-7(7 ■ 



13 



(7-h 2cos^6i- 1) 



(1)|2 



/3 -f 7(7 — /32 



2UJ (/3 + 7(7)(l + fc2)2 



-K/cfg [p{l + k^ + 2saP) +(7 7(7-l)(l + /c2)2i 



_ f20)|,(l)|2 



= 7 



\4>\ 



(1)|2 _ ^(20) I ,(1)|2 



(s/3B -H2ct) 



(1)|2 _ (20) I ,(1)|2 



and 



'v,B~2^-^±^cose 



I0i 



(1)|2 



(25) 



(26) 



It is expected, and indeed verified by a tedious yet straightforward calculation, that upon setting cr = 0, s = — 1 in 
expressions (|24|) and (|25|l . one recovers exactly Eqs. (44) - (49) in Ref. [given (42) therein]. 



Notice, for rigor, that for 'vanishing obliqueness' i.e. if 6* ^ 0, one obviously has k • u| 



(") 



k (by definition) , 



implying the condition: C2 



nl) 



kc"^''^ (for 9 —!■ 0) which is indeed satisfied for all n, l,hy the above formulae. 
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E. Derivation of the Nonlinear Schrodinger Equation 



Proceeding to the third order in e (n — 3), the equation for / = 1 yields an exphcit compatibihty condition to be 
imposed on the right-hand side of the evolution equations which, given the expressions derived previously, can be cast 
into the form 



(27) 



where = '^i^'' denotes the amplitude of the first-order electric potential perturbation; coefficients Ai^2.3 are to be 
defined. Now, multiplying by iA^^, we obtain the familiar form of the Nonlinear Schrodinger Equation 



9-0 



0. 



(28) 



Recall that the 'slow' variables {C,'''} were defined in 0. 

The dispersion coefficient P = —A2/A1 is related to the curvature of the dispersion curve as P 



1 d^uj 

2 alj 



w"(fc) cos^ei +u'{k) 



sin^ 9 



the exact form of P reads 



1 1 

]32^ 



vi - (i/i +3^^^) cos^6' 



(29) 



where we have defined 



,3 3 /3-H7(T(3-P)(l-ffc^) 
3[/3 + 7a(l + fc2)]4 



(30) 



Note that, just like vi defined above, 1^2 ^ 1 when cr — > 0; see that relation (51) in Ref. 20j is recovered from H29|) 
in this case. If, furthermore, we set [3 — 1 {m addition to cr = 0) in all expressions describing our dispersion law i.e. 
(P|l . (EU, lES|) above, we obtain respectively (3), (11), (4) in Ref. l3- 

It seems appropriate, here, to point out the qualitative difference between P given in (|29|l as compared to relevant 
previous expressions: the existence of a may affect the sign of the P coefficient. For instance, taking cr = (i.e. 
ui = V2 = P '^s readily seen to be negative for parallel modulation, i.e. setting = 0; however, for cr ^ this is 
no longer the case, since P changes sign at some critical value of k (to see this, study the sign of V2 versus k ii])- 
Furthermore, a similar remark holds for the effect of an oblique modulation on the sign of P; we will come back to 
this subtle point in the next subsection. 

The nonlinearity coefficient Q = —A^/Ai is due to the carrier wave self-interaction. Distinguishing different 
contributions, Q can be split into five distinct parts, viz. 



= Qo + Qi + Q2 + Q.3 + Qa 



(31) 



reflecting the similar structure of A3 



Az = A3 



(0) 



A 



(1) 



A 



(2) 



A 



(3) 



a: 



(4) 



(32) 



In order to trace the influence of the various parameters, let us define all quantities in full detail. First, a'"^^ (as well 



as Qo 



-Ag^^/Ai) is related to the self-interaction due to the zeroth harmonic, viz. 



a: 



(0) 



-I3k'{c\ 



(11) J20) 



+ C. 



(11) J20) 



^ ~ sujlak C4 c 



2^(11)^(20) 



UJ (1 



k^ 



J11)J20) 



(33) 



while Ag^-* (related to Q2 = —A^^^'/Ai) is the analogue quantity due to the second harmonic 



- _4(2) 



A^f=-Pk'{c^^ 
„("0 



„(11)J22) 



s ^ 2a fc2 cf^cf ^ -c.(l + fc2) 411)422) 



(34) 



All coefficients Cj were defined previously. Now, Qi 
cubic term in (|HJi) (often omitted in the past) 



-A^^'' /Ai is simply the nonlinearity contribution from the 



A 



(1) 



(35) 
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Finally, ylg^-* (related to Q3 = —A'f^/Ai) is the {a- related) result of the third line in 



l(3) 



l(3) 



while Ag'*'' (and Q4 = — Ag'^-'/Ai) is due to the last two lines in 



We note that Ai is everywhere defined as 
i.e. by using (fTS|l 



A, 



1 1 /uj' 



(36) 



(37) 



(38) 



(39) 



(reducing to: A^^ ^ Tp W ^'^^ cr = 0). Remember that Q3 and Q4 are plainly absent from the previous results in 
Ref. 20] (i.e. for a — 0) and so is, in fact, Qi. 

Substituting from the expressions derived above for the coefficients cf^^^ and re-arranging, we obtain 



1 



1 



2UJ /?2 (l + /c2)2 p^^fj^y2 



(3k^ 



f3[3 + 6k^ + 4k'^ + k^ + 2af3{s {2k^ + 3) + 2 a w^)] 
+ 7(7 [ (7 + 1) (1 + k^ f + 2a(3{~2a/3 + S7 (1 + k^)] 
+ [(3 {2 + 4/j2 + 3/c4 + + 2saP) + 2ja{l + k^ {1 + k^ + saP)] cos 29 
+ 2{l + kY{P + lcr)(^^ cos^O 



n = 



+ k{i + k^) 



/3/fc2 + (1 + A:2) 



— X 



P{l + k^ + 2saf3) + 7 (7 - 1) a (1 + k^ 



cos ( 



2UJ (l + fc2)2' 
1 1 1 



(40) 
(41) 



12 /?3 LU fc2(l + fc2)2 

2(3 k^ 5saP^{l + k^f + 2a^P^ + 2 7 V (1 + fc^)^ (1 + 4fc2) 
+ /3 (1 + k^f (3 + 9fc2 + 2 s a 7^ a) 



+ (1 + P)3^2 



/5 (3 + 9fc2 + 6fc4 + 2saP) + 2 7^ cr (1 + fc2)2 (1 + Ak^) 



(42) 



Finally, the coefficients Qz = —A^^^/Ai and Q4 = — ^g^-*/^! can be directly computed from 13t)|) - (|38|) above; the 
lengthy final expressions are omitted here. 

Once substituted in H31|l . these expressions provide the final expression for the nonlinearity coefficient Q. One may 
readily check, yet after a tedious calculation, that expressions (jlU)! and reduce to (53) and (54) in Ref. for 
a — Q. However, the remaining coefficients Qi, Qs, Q4 were absent in all previous studies of the DA waves, to the 
best of our knowledge. Their importance will be discussed in the following. Note that Qi, Q2 do not depend on the 
angle 9. 
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F. Behaviour of coefficients for small k 

A preliminary result regarding the behaviour (and the sign) of the NLSE coefRcients P and Q, at least for long 
wavelengths, may be obtained by considering the limit of small fc ^ 1 in the above formulae. 

The parallel {9 = 0) and oblique {9 ^ 0) modulation cases have to be distinguished straightaway. For small values 
of fc (fc <C 1), P is negative and varies as 

in the parallel modulation case (i.e. 9 — G), thus tending to zero for vanishing fc, while for 7^ 0, P is positive and 
goes to infinity as 



for vanishing fc. Therefore, the slightest deviation by 9 of the amplitude variation direction with respect to the 
wave propagation direction results in a change in sign of the dispersion coefficient P. Given the importance of the 
coefficient product PQ (to be discussed in the next Section) , one may wonder whether this is sufficient for the stability 
characteristics of the DA wave to change. Let us see what happens with the Q in the limit of small fc. 

For all cases, Q varies as ~ 1/fc for small fc <C 1 js^l; the exact expression in fact depends on the angle B. In the 
general case (0 7^ 0), the result reads 

QL^o ^ ^^'"^ + + ^^'""^ ^^'"'^ + 3) + 7 (7 + 1) ^] ^ • (45) 

A careful study shows that Q is negative, in fact, for all possible values of the physical parameters of interest (i.e. a, 
/?, 7, (T - all positive - and s ± 1). For vanishing 9, however, the approximate expression for Q, yet apparently quite 
similar, is now positive, i.e. 



9=0 



^ ^ [/3(2sa/3 + 3) + 27cr][/3(2sa^ + 3)+7(7 + l)CT]i. (46) 



In conclusion, both coefficients P and Q change sign when 'switching on' theta. Indeed, obliqueness in modulation 
is expected to inffucncc the stability profile of the system; this point seems to confirm (and complete) the general 
qualitative arguments put forward in Ref. [T^ for the ion acoustic wave in an electron ion plasma without dust. 
Nevertheless, at all cases, the product of P and Q is negative for small k, ensuring, as we shall see in the following 
section, stability for long perturbation wavelengths. As a by-product of this analysis, we see that taking into account 
Qi, Q3 and Q4 does not seem to influence the dynamics in the low wavenumber k parameter range. 

IV. STABILITY ANALYSIS 

The standard stability analysis [8, 34] consists in linearizing around the monochromatic (Stokes's wave) solution of 
the NLSE ^ 

(notice the amplitude dependence of the frequency) by setting 

ip = il^o + e-ipi , 

and taking the perturbation ■01 to be of the form: -01 — ^1,0 e*'-*'^^'^^'' + c.c, (the perturbation wavenumber k and 
the frequency uj should be distinguished from their carrier wave homologue quantities, denoted by k and uj). Now, 
substituting into l|28|) . one readily obtains the nonlinear dispersion relation 

co' = P'P(P -2^\i,,A. (47) 



One immediately sees that the wave will be stable for all values of k if the product PQ is negative. However, for 

^1 



positive PQ > 0, instability sets in for wavenumbers below a critical value kcr = \/2^\ipi^o\, i.e. for wavelengths above 
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In brief, we see that the instabihty condition depends only on the sign of the product PQ, which can now be studied 
numericaUy, relying on the exact expressions derived in the preceding Section. 

In the contour plots presented below (see figures El), we have depicted the PQ = boundary curve against the 
normalized wavenumber k/ku (in abscissa) and angle 9 (between and tt); the area in black (white) represents the 
region in the (k — 9) plane where the product is negative (positive); instability therefore occurs for values inside the 
white area. We have considered values of the wavenumber k between zero and upto 4 times the Debye wavenumber 
ko (yet mostly focusing our attention on the low k region). Pitch angle 9 is allowed to vary between zero and 7r/2; 
as a matter of fact, all plots are ^- periodic, i.e. symmetric upon reflection with respect to either the 6* = or the 
9 = ^ lines. We have chosen a fixed set of representative values: a = 5 ■ 10~^, a' = 2a^/3 « 1.6 • 10^^ and f3 « 100, 
corresponding to Zd/Zi = 10'^ and T^/Ti — 10 (we have taken 7 = 2, cr = 1 for the plots). 

For negative dust (s — ~1; see fig. ^ the product possesses positive values for angle values between zero and 
9 K, 51°; we see that instability sets in above a wavenumber threshold which is clearly seen to decrease as the 
modulation pitch angle 9 increases from zero to approximately 17 degrees, and then increases again up to « 51°. 
Nevertheless, beyond that value (and up to 7r/2) the wave remains stable; this is even true for the wavenumber regions 
where the wave would be unstable to a parallel modulation: see e.g. the interval where 9 = Q and k/kjj £ [1.0,3.6] 
approximately, in figure ^ The inverse effect is also present: even though certain k values correspond to stability for 

= 0, the same modes may become unstable when subject to an oblique modulation {9 ^ Q); this is mostly true for 
long wavelengths (small k). Notice the periodicity with respect to 9. 

A similar behaviour is witnessed in the case of positive dust (s = +1; see fig. 2), yet the instability threshold kcr 
for a given value of 9 is quite higher: positive dust rather appears to favour stability. 

In all cases, the wave appears to be globally stable for large angle 9 modulation (between 0.9 and 7r/2 radians, 
i.e. 51° to 90°) and unstable for smaller values of 9. For parallel modulation (9 = 0), the sign of the product PQ is 
basically opposite to that of Q, since P < for all values of fc; the wave is then stable for large wavelengths A ^ Ad 
(i.e. for k/ko <^ 1), and potentially unstable for higher values of k (a similar qualitative behaviour has been reported 
for the ion-acoustic wave case (i.e. without dust) [11 - 15]. 

A final word is in row, concerning the effect of taking into account the pressure evolution equation Q , often omitted 
for simplicity. Given the above results, this amounts to wondering what the difference would be, should we simply set 
cr = in expression for P and thus plainly omit Q3 and Q4, defined above. A qualitative answer is attempted 
in figure 13 where we have depicted the PQ product in this case. The qualitative results obtained so far do not seem 
to be strongly modified, at least not for low values of k (say, below k « 2kD) and definitely not as far as the angle 
dependence of stability is concerned. The difference in stability regions obtained for higher k is rather negligible for 
long wavelengths (say, below k « l.^kn)- Nevertheless, including the pressure equation in the description seems to 
describe the problem in a more precise manner, and also somewhat restricts the instability region, since stability is 
now predicted for short wavelengths (above, say, k « 3.6fc£)) and low 9; compare figs. ^,|2K to|5^,|ntj respectively. 



Let us discuss the possibility of the existence of localized excitations in our system. The NLSE (|28|l is known to 
possess distinct types of localized constant profile (solitary wave) solutions, depending on the sign of the product PQ. 
We shall now briefly outline the method employed to derive their form and discuss their relevance to our problem. 
Following Ref . [33 , we may seek a solution of Eq. H28|) in the form 




where p, a are real variables which are determined by substituting into the NLSE and separating real and imaginary 
parts. The different types of solution thus obtained are clearly summarized in the following paragraphs. 



V. NONLINEAR EXCITATIONS 




(49) 



A. Bright solitons 



For PQ > we flnd the (bright) envelope soliton |37| 





(50) 
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representing a localized pulse travelling at a speed u and oscillating at a frequency Vt (for u = 0). The pulse width L 
depends on the (constant) maximum amplitude square po as 

\ Q Pa 



B. Dark solitons 



For PQ < we have the dark envelope soliton [hole) |31 



P = Pi 



1 — sech 



21 C - 



= pi tanh 



21 C-UT 



(52) 



representing a localized region of negative wave density (shock) travelling at a speed u. Again, the pulse width 
depends on the maximum amplitude square pi via 



L'=J2 



P 

Qpi 



(53) 



C. Grey solitons 

It has been shown in Ref. [s^ that looking for velocity-dependent amplitude solutions, for PQ < 0, one obtains 
the grey envelope solitary wave 



P 

e 



P2 

1 

2P 



1 2 1,21 C-ur 
1 — a seen 



VoC 



- 2PQp2 ) T + eio 



a tanh( '' 



seen 



1/2 ' 



(54) 



which also represents a localized region of negative wave density; Oio is a constant phase; S denotes the product 
S — sign P x sign[u— Vb). In comparison to the dark soliton (|52|l . note that apart from the maximum amplitude 
P2, which is now finite (i.e. non-zero) everywhere, the pulse width of this grey-type excitation 



i"- W2 



P 



QP2 



now also depends on a, given by 



1 



2Pg P2 



in' Vi) < 1 



(55) 



(56) 



(PQ < 0), an independent parameter representing the modulation depth (0 < a < 1). Vq is an independent real 
constant which satisfies the condition [3al 



Vo - ^2\PQ\p2 <u<Vo + ^j2\PQ\p2 ; 

for Vq = u, we have a = 1 and thus recover the dark soliton presented in the previous paragraph. 

Summarizing, we see that the regions depicted in figs. |21 13 El actually distinguish the regions where different 
types of localized solutions may exist: bright (dark or grey) solitons will occur in white (black) regions (the different 
types of NLS excitations are exhaustively reviewed in [SS])- Furthermore, soliton characteristics will depend on the 
dispersion laws via the P and Q coefficients; for instance, regions with higher values of P (or lower values of Q) - see 
figs. 1301- will support wider (spatially more extended) localized excitations. 
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VI. CONCLUSIONS 

This work has been devoted to the study of the conditions for occurrence of the modulational instabihty of the 
dust-acoustic waves propagating in an unmagnetized dusty plasma. Considering the Poisson-moment equations for 
the dust and aUowing for modulation to occur in an oblique manner, we have shown that the DA wave modulational 
instability depends strongly on the angle between the propagation and modulation directions. As a matter of fact, 
the region of parameter values where instability occurs is rather extended for angle 9 values up to a certain threshold, 
and, on the contrary, smeared out for higher 9 values (and up to 90 degrees, then going on in a ^ - periodic fashion). 

Furthermore, we have studied the possibility of the formation of localized structures (solitary waves) in the system. 
Distinct types of localized excitations (envelope solitons) have been shown to exist. Their type and propagation 
characteristics depend on the carrier wave wavenumber k and the modulation angle 9. 

Summarizing our results, we have seen that 

(i) obliqueness in the amplitude modulation direction has a strong influence on the conditions for the modulational 
instability to occur: regions which are stable to a parallel modulation may become unstable when subject to an 
oblique modulation, and vice versa; 

(ii) large-angle modulation seems to have a stabilizing effect; on the contrary, small-to-medium angle (say below 50 
degrees) modulation enhances instability; 

(iii) DAW-related localized excitations may appear and propagate in a dusty plasma; modulationally stable (unstable) 
(fc, 9) regions support envelope solitary waves of the bright (dark) type; 

(iv) the type and characteristics of the latter (localized modes) depend on the value of 9: for given low k, dark solitons 
(or holes) are wider as 9 becomes higher (see fig. 0J|; for higher fc, bright (dark) solitons become narrower (wider) as 
9 increases; finally, for given 9 values below (above) a threshold of, say, 51 degrees, bright (dark) excitations will be 
narrower (wider) for higher k (see fig. 0}; 

(v) comparing the positive (s = -1-1) to negative (s = —1) dust cases, we have shown that positive dust enhances 
stability and rather favours dark-type excitations (hole solitons); furthermore, low k dark envelope solitons appear to 
be narrower with positive dust; for higher k there is practically no qualitative difference between the two dust charge 
sign cases; 

(vi) As a final comment, let us point out that taking the dust pressure equation ^ into account, we have obtained a 
wider stability region for small 9 values, yet only for high wavenumbers. The existence of dark-type localized envelopes 
of high k modes subject to slightly oblique (low 9) modulation is thus predicted; cf. figs. to|Sli,|ntL, respectively. 
However, for wavenumbers below, say, k = kjj, there is no qualitative difference due to the consideration of |(2J). 

Our aim has been to put forward a model study of the DAW modulation which is generic, i.e. incorporating 
several previous descriptions, which may be recovered for different choices of the physical parameters involved in 
the formulation. Dust charge was assumed to be constant and the plasma geometry was taken to be Cartesian and 
infinite, for simplicity. Thus, our work complements the investigation by Tang and Xue '38*1 who examined only the 
modulational instability of DAWs against oblique modulations, including an ad hoc charging equation and a specific 
form of the adiabatic law for warm charged dust grains which are negatively charged. The present paper, on the 
other hand, discusses the multi- dimensional modulational instabilities of dust acoustic waves in plasmas containing 
both negatively and positively charged dust grains, as well as provides a detailed discussion of various types of dust 
acoustic envelope solitons and their respective parameter regions of existence, leaving the choice of the value of the 
parameter 7 = Cp/cy free in the algebra. 
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Figure captions 

Figure 1: 

(a) The coefficient product PQ = curve is represented against normalized wavenumber k/ku (in abscissa) and 
angle 9 (between and tt); the area in black (white) represents the region in the (fc — 9) plane where the product 
is negative (positive); instability therefore occurs for values inside the white area. This plot refers to negative dust 
charge (s = — 1). (b) A close-up plot near the origin. 

Figure 2: 

Same as in figure 1, for positive dust charge (s = +1). Notice that the stability region close to the origin gets 
narrower. Positive dust charge seems to favour stability. 

Figure 3: 

(a) The curves for constant values (contours) of the dispersion coefficient P are represented against normalized 
wavenumber k/ko (in abscissa) and angle 9 (between and tt/2); In ascending order (from bottom to top), the curves 
correspond to P = —0.4, —0.3, 0.3, 0.4; P clearly increases with 9, for a given wavenumber k. The parameters used 
for this plot are as defined in fig. 1. (b) A similar contour plot for the nonlinearity coefficient Q. In descending order 
(from top to bottom), the curves correspond to Q = -0.003, -0.0025, -0.002, -0.002, 0, 0.001, 0.002; Q decreases 
with increasing 9, in this region. Remember that (the part of) these curves falling inside the instability region (white 
sector in fig. ^ is related to the instability growth rate a via H48(l . Values of fT/|'0i,oP above a certain value are to be 
excluded, since they would fall inside the stability (black) region: this element is absent in fig. ^ This plot refers to 
negative dust charge (s = — 1). (c) The analogous contour plot (same values as in (b)) for the nonlinearity coefficient 
Q in the positive dust charge (s = +1) case; Q takes higher values here (cf. (b)), for given {k,9), leading to a more 
extended stability region for large wavelengths (fc <C ko). 

Figure 4: 

Contours of the ratio P/Q (whose absolute value is related to the soliton width; see H51|l . (|53|l ') are represented 
against normalized wavenumber k/k^ (in abscissa) and angle 9 (between and 7r/2); In descending order, starting 
from above, the curves correspond to P/Q = —20, —10, —5, —1, 0, 1, 5, 10, 20; the value of P/Q decreases with 9, 
for a given wavenumber k above ko, so higher 9 seems to favour narrower (wider) bright- (dark-) type excitations. 
The same qualitative behaviour was obtained for positive dust charge i.e. s = +1 (not depicted, for the difference 
was unimportant). 

Figure 5 

The product PQ, as in fig. QJl, as results from the pressure equation Q being omitted. Comparing to fig. ^ notice 
that there is practically no qualitative difference for low k and for high 9; however, predicted behaviour changes above, 
say, k Ri 1.5 ku. This plot refers to negative dust charge (s = —1). 

Figure 6 

Similar to O but for s = +1 (positive dust charge) i.e. as in fig. 12 but omitting the pressure equation Once 
more, the change in the qualitative analysis does not appear to be dramatic. 
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